Suppression of non-adiabatic phases by a non-Markovian environment: 
easier observation of Berry phases. 
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We consider a two-level system coupled to a highly non-Markovian environment when the coupling 
axis rotates with time. The environment may be quantum (for example a bosonic bath or a spin 
bath) or classical (such as classical noise). We show that an Anderson orthogonality catastrophe 
suppresses transitions, so that the system's instantaneous eigenstates (parallel and anti-parallel to 
the coupling axis) can adiabatically follow the rotation. These states thereby acquire Berry phases; 
geometric phases given by the area enclosed by the coupling axis. Unlike in earlier proposals 
for environment- induced Berry phases, here there is little decoherence, so one does not need a 
decoherence-free subspace. Indeed we show that this Berry phase should be much easier to observe 
than a conventional one, because it is not masked by either the dynamic phase or the leading non- 
adiabatic phase. The effects that we discuss should be observable in any qubit device where one 
can drive three parameters in the Hamiltonian with strong man-made noise. 

PACS numbers: 03.65. Vf, 03.65.Yz, 85.25.Cp 



I. INTRODUCTION 

Noise is typically a huge inconvenience when trying to 
control the state of a quantum system, since it causes dis- 
sipation and decoherence. However one can ask if noise 
could be used to coherently control the quantum sys- 
tem in a manner that cannot be achieve using traditional 
Hamiltonian manipulation. 

In this context we consider the Berry phase [l| ; a geo- 
metric phase associated with adiabatic evolution. When 
one rotates the parameters in a system's Hamiltonian 
around a closed loop slowly enough that an cigenstate 
adiabatically follows the change, then that state acquires 
a Berry phase, "J'bp- This phase depends on the geom- 
etry of the parameters' path but not on how that path 
is followed. For a spin-half in a slowly rotated magnetic 
field, $bp = — s z A, where A is the solid-angle enclosed 
by the field and s z is the quantum number of the spin 
along the field axis. Berry phases occur in many quan- 
tum systems @, H, H|, they were recently observed in 
superconducting qubits and their noise-dependence 
has been investigated using cold- neutrons [|| . They have 
potential applications in both quantum computation [?J 
and metrology Q, because it is argued that one of the 
most accurate ways to change the phase of a state is to 
rotate the Hamiltonian's parameters round a loop which 
encloses a given solid angle. 

However there is a practical problem with observing 
the Berry phase; one does not observe it "alone" . If the 
parameters of a system's Hamiltonian complete a closed 
loop in a time i p , then the total phase acquired by an 
eigenstate is [§) 

$total = $dyn + *BP + + $ NA + " " " > C 1 ) 

where the dynamic phase $dyn oc Et v , the Berry phase 
$bp oc (Et p )°, and the non-adiabatic (NA) correction 
^na ^ (Etp) -1 *! w ith E being a system energy scale fl0| 
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FIG. 1: Sketches of the two situations we consider in this 
work, (a) A Berry phase is created in the two-level system 
by coupling it to a quantum environment. The environment 
coupling is along the axis e(t), which is slowly rotated around 
a closed loop, (b) A Berry phase is created in the two-level 
system, by coupling it to classical non-Markovian noise, A(t). 
Here the noise is along the slowly rotating e(t)-axis. In both 
(a) and (b), the Berry phase is given by the solid-angle en- 
closed by e(t). 



(such as the gap to excitations). As the second term in 
the t~ ^expansion of $totai, the Berry phase is difficult 
to isolate, and thus hard to utilize for quantum compu- 
tation or metrology. One must make t p large to suppress 
the non-adiabatic terms. This makes 'E'dyn 3> 3>bp, so 
one must subtract off $dyn (using a spin-echo trick Q 
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or degenerate states [ll|) with extreme accuracy. Con- 
cretely, if one wants <I>bp to an accuracy of one in 10 3 , 
one requires Et p ~ 10 3 . Thus $dyn ~ 10 3 $bp, so one 
must subtract off $dyn with an accuracy of one in 10 6 . 
Furthermore, as t p must be long, any device is slow to 
operate and leaves lots of time for decoherence to destroy 
the phase. 

In this work, we analyze a system coupled (via a term 
Hint in the Hamiltonian) to a highly non-Markovian envi- 
ronment with N modes (see Fig.QJ,). Such environments 
are known to strongly renormalize the system dynamics; 
inducing a type of Anderson orthogonality catastrophe 
fl2{ , as shown by Leggett et al usin g a method that they 
called adiabatic renormalization [13|,llJ]. For large N, we 
show that when 7ii n t is changed slowly the instantaneous 
cigenstatcs of Hint adiabatically follow the change. As- 
suming the system has no intrinsic dynamics, the total 
phase acquired by an instantaneous eigenstate when the 
parameters of Ti.- lnt are rotated around a closed loop is 

*total = 3>BP + *n! + ' ' ' • (2) 

Both the dynamic phase, ^dyn, and the first non- 
adiabatic phase, $na> are absent. As usual the Berry 
phase, $bp, is half the solid-angle enclosed by the pa- 
rameters of the Hamiltonian (in this case the solid-angle 
enclosed by the coupling axis in 7ii n t)- 

Crucially, it is much easier to accurately observe the 
Berry phase in a situation where $totai is given by Eq. ([2]) 
in place of Eq. ([T]). There is no dynamic phase to sub- 
tract off, and the leading non-adiabatic phase is ab- 
sent. If we again assume we want to get the Berry 
phase with a accuracy of one in 10 3 , then we require 
that = {E2t p )~ 2 ~ 10~ 3 , which simply means that 

E 2 t p ~ 10 3 / 2 <~ 31 (where E 2 is the energy scale in $^ 2 |). 
This is a vastly less difficult to achieve than the condi- 
tions discussed above for a conventional Berry phase. 

We note that there is another geometric phase - 
the Aharonov-Anandan phase [l5| — which has no non- 
adiabatic corrections. Thus one may expect that one can 
use it to avoid the above problems with non-adiabatic 
phases. However the Aharonov-Anandan phase has dif- 
ferent properties from the Berry phase, which make it 
more difficult to calculate in many experimental situ- 
ations. We reserve discussion of this difficulty to Ap- 
pendix 13 and here consider only the Berry phase. 

The effects that we discuss above for a quantum en- 
vironment, are equally relevent for a system coupled to 
classical noise. Classical noise is known to be equivalent 
to a hig h-temperature environment of harmonic oscilla- 
tors [16| ■ We use this equivalence to show that driving a 
system with highly non-Markovian noise (such as strong 
high-frequency noise) whose axis is slowly rotated (see 
Fig. [TJd) causes instantaneous eigenstates to acquire a 
phase given by Eq. ([2]). 

For cither model (quantum environment or classical 
noise) , there will be some dephasing of the system state, 
however it originates from the non-adiabatic rotation of 



the coupling axis. Thus if the rotation is performed in 
a time t p , the dephasing goes like exp[— (£^ii p ) -1 ] so the 
longer the experiment the less the dephasing. For all 
t p > E± , the dephasing is weak-enough to clearly mea- 
sure the phase. Thus there is no need to work in a 
decoherence-free subspace to observe this environment- 
induced Berry phase (unlike that in Ref. [2^| which we 
discuss in Section II Dl below) . The only difference be- 
tween a low-temperature quantum environment and clas- 
sical noise is in this dephasing. For low-temperature 
quantum environments, E\ grows exponentially with the 
strength of the environment coupling, while for classical 
noise, E\ goes like the squareroot of the noise power. 
Thus for strong coupling to a quantum environment one 
can easily be in a situation where the (-Eii p ) _1 is irrele- 
vant, then it turns out that the leading contribution to 
dephasing goes like (E^tp)^ 3 (where we believe that E$ 
is not exponentially suppressed) , so it scales to zero much 
faster with increasing t p than for classical noise. How- 
ever we emphasize that the dephasing is already weak 
for classical-noise, so this fact that it is much weaker for 
a low-temperature quantum environment is not central 
message of the work we present here. For all practical 
purposes, in the context of the problems we consider 
here, classical noise and quantum environments induce 
the same effects. 



A. Two-level system with non-Markovian quantum 
environment 

We consider a spin-half coupled to an environment via 
a coupling axis, e(f), which is rotated around a closed 
loop (see Fig. QJt) . The total Hamiltonian for the system 
(sys) and its environment (env) is 

}), (3a) 

H int (t) = -\a ■ e(t) K M + (3b) 
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with & = (a x , a y ,a z ) being the vector of Pauli matrices. 
The unit vector e(t) is slowly rotated around a closed 
loop. The operators and cij, create and destroy the 

jth excitation of the environment, where Ti env ({a^, cij}) 
is such that the jth excitation has energy Cij . Note that 
this is a model of a degenerate two-level system coupled 
to an environment, because we require that the two-level 
system's Hamiltonian is zero if one takes 7ii n t — > (we 
relax this requirement in Section IVI[) . 

We assume that the environment has a spectrum such 
that it is highly non-Markovian; by which we mean 
that the environment has a significant effect on the sys- 
tem's dynamics on a timescale much less than the en- 
vironment's memory time, r mom . This memory time is 
defined as the timescale on which X^("]( r )%(0)) de- 
cays, where aj(r), a,j(r) are creation and annihilation 
operators in the interaction picture, so that a At) = 
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cxp[-i7Y cnv r]a] cxp[i7Y cnv r]. 

Here we say that any zero temperature environ- 
ment is highly non-Markovian when the dimensionless 
environment-coupling parameter 



G dis = f! m 2 K ] » !. 



(4) 



where f2 m is the average frequency of the environment; 
defined by Q m = KjQjj / fej K f) ■ For an envi- 
ronment of harmonic oscillators this can be generalized 
to arbitrary environment temperature, T; such an envi- 
ronment is highly non-Markovian if 



G dis = ft m 2 ^!^coth[/3%/2] » h 



(5) 



where /3 = (k^T)" 1 . We have in mind an environment 
containing so many modes with different frequencies that 
the sum can be written as an integral as in Eq. f|23|) . 

To get a feeling for such an environment, it is worth 
considering an environment at zero temperature with N 
modes each coupled to the system with a strength K 7 
and with frequencies spread around an average of Q m . 
For simplicity we assume that the spread of frequencies 
is also of order fi m , in such a case Gdis — NK 2 /Q 2 n . 
One can then see that Gdis "C 1 corresponds to an ap- 
proximately Markovian environment by looking at the 
derivation of the Bloch-Rcdfield master equation 



which has a Markovian form (see for example Ref. l2lf h 
There the memory time r mcm ~ ^i^ 1 , while the rate at 
which the environment affects the system's dynamics is 
r ~ NK 2 r mcm . The Bloch-Rcdficld approximation is 
then applicable when rr mem <S 1; then spin relaxation 
and dephasing times (known as T\ and Ti respectively) 
are of order T _1 so r mcm <C Ti.2- This corresponds to 
Gdis < 1. In contrast the Lindblad equation master 
equation [1?], [H, H3 is only applicable when r mcm — > 
(strictly Markovian environment) with NK 2 scaled so 
that r remains finite. This therefore corresponds to 
Gdis = 0. 

In this work we consider exactly the opposite limit, 
Gdis 3> 1. One might guess that this corresponds to 
dephasing so strong that it occurs in a time shorter than 
the environment memory time. However in this regime 
relaxation and dephasing rates are not given by the above 
r. Instead there is extremely strong renormalization of 
the system's dynamics with relatively weak relaxation 
and dephasing. 



B. Environments not made of harmonic oscillators 

Most works on systems coupled to environments have 
considered that the environment is made up of har- 
monic oscillators. We do the same here, treating Ti. cnv = 
QjCLjCLj (up to an irrelevant constant that we neglect). 



This enables us get closed form expression for finite tem- 
peratures without too much difficulty. However, we em- 
phasis that the low temperature limit of our results are 
applicable to any large environment, for example an en- 
vironment made of two-level systems. 

To see why this is so, consider any large environ- 
ment which is initially in its ground-state. The analy- 
sis that follows in this article assumes that the system 
only weakly interacts with each degree-of-frccdom in the 
environment, The effect of the environment on the sys- 
tem is none the less strong because there are so many 
degrees-of- freedom in the environment. Thus the inter- 
action with each environment degree-of-freedom can be 
treated to lowest order. This means that the interaction 
can only take environment mode j from its ground state 
to its first excited state. In this case the only property 
of environment mode j that is relevant to the system's 
evolution is the gap between its ground state and its first 
excited state, which we define as fij. Thus, for temper- 
atures such that e _/3Sl3 <C 1 for the vast majority of j, 
the nature of environment modes is strictly irrelevant,; 
they could be harmonic oscillators (with an infinite lad- 
der of equally spaced excited states), spin- halves (with 
no states beyond the first excited state), or any other 
quantum systems. 

Thus one can expect that the /3r2 m ^ I limit of all re- 
sults in this article will apply to arbitrary environments, 
while the results for arbitrary (3 only apply to an envi- 
ronment of harmonic oscillators. 



C. Two- level system with non-Markovian classical 
noise 

One can map the dynamics of a system coupled to an 
environment of quantum oscillators at infinite tempera- 
tures onto the ensemble-averaged dynamics of a system 
coupled to classical noise [l6[. Thus we also analyze the 
problem of a spin subject to classical noise (with TV com- 
ponents) when the noise axis is rotated (see Fig. [Ud). 
In Section [V] we discuss what happens to a spin whose 
Hamiltonian is given by 

H c i = -±a-e(t)J2 A j cosiftjt + fa), (6) 
i 

when the unit vector e(t) is slowly rotated around a 
closed loop, and Aj,4>j are random. The amplitude Aj 
is gaussianly distributed with variance (A 2 ), and <pj is 
uniformly distributed over all angles from to 2ir. 

The mapping from the quantum environment to clas- 
sical noise tells us that the average system dynamics un- 
der 7i c i is given by taking the infinite temperature limit 
(/3 = (/sbT) -1 — ► 0) of the system's dynamics under 
Wsys&env when the environment consists of harmonic os- 
cillators and K 2 is replaced by (3£lj(A 2 ). Thus for such 
classical noise we can define the noise-coupling parameter 



Gdi, 



o- 



(7) 
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Here, we consider only highly non-Mar kovian classical 
noise for which Gdis ^ 1- 



D. Prior works Berry phases with environments 
and non-adiabatic corrections. 



When one performs a Berry phase experiment on a 
system (rotating its Hamiltonian around a closed path), 
one can rarely ignore the coupling to an environment (or 
noise). Despite hints to the contrary [22], noise-induced 
fluctuations of $dyn typically lead to a dcphasing time, 
T2, on which all pha se information (including the Berry 
phase) is lost [23l. |24|. [25j (here we do not consider the in- 
triguing use of repeatable noise in Ref. Q ) . Thus a Berry 
phase can only be observed if one is able to adiabaticaUy 
rotate the Hamiltonian around a closed loop in a time 

;$ T 2 . Usually this requires T 2 > t p ^> E~ x , where E 
is the energy gap to system excitations frol | . The environ- 
ment also modifies the Berry phase [23| ; the modification 
is geometric (quadrupole-like) and complex [26|, [27j (see 
also [28j|). with its imaginary part being geometric de- 
phasing. 

Ref. [29[ made the remarkable observation that a time- 
dependence in the coupling to an environment could 
also generate a Berry phase (see also [30|, HU). Dissipa- 
tive processes (at a rate T 2 _1 ) cause the system state to 
adiabaticaUy follow the time-dependence of the system- 
environment coupling, thereby acquiring a Berry phase. 
In Ref 29], unlike in our Eq. ([2]), there are non-adiabatic 

corrections of the form ~ {tp/T2)~^ for all integer 
fi > 1. To avoid the non-adiabatic effects, one requires 
that t p 3> T2, which means that dephasing occurs long 
before the rotation is completed. Thus this Berry phase 
can only be observed for states in a decoherence-free sub- 
space [29l |. 

All the works listed above considered Berry phases in 
systems coupled to environments or classical noise that 
were strictly or approximately Markovian. In contrast, in 
this article we consider a highly non-Markovian environ- 
ment, which strongly renormalizes the system dynamics 
without causing significant dissipation. The adiabatic 
evolution is ensured by this renormalization (not dissipa- 
tion), and it leads to Eq. ([2]). The relative lack of dissipa- 
tion means that there is no need to use a decoherence-free 
subspace. 

Neither the Lindblad nor Bloch-Redfield methods, 
used to study Berry phases in Refs. [29|,[3(| and Ref. [3lJ 
respectively, can capture the strong renormalization that 
occurs for Gdis ^ 1- The Lindblad master equation 
applies for Gdis = 0, while the Bloch-Redfield master 
method applies for Gdis <C 1, which corresponds to only 
a very small linear renormalization effect. 

We conclude by mentioning a number of works in 
chemical physics. Berry phases occur in molecular dy- 
namics because there is a separation of timcscales; the 
nuclei are heavy and move slowly, while the electrons 
are light and move fast. This makes it natural to per- 



form a Born-Oppcnheimer decoupling of the fast and 
slow degrees-of-freedom. It is well-known that such a 
de-coupling can lead to a Berry phase, see for example 
Refs. @, There are also non-adiabatic corrections to 
this Berry phase which come from violations of the Born- 
Oppcnhcimer decoupling. These have been well studied; 
for an older review see Ref. 1321] and references therein, 
for more recent work see Ref. [33| . However all the works 
that we are aware of, consider only a relatively small 
number of degrees-of-freedom (for example a tri-atomic 
molecule, containing three slow nuclear degrees of free- 
dom and up to three fast electron degrees-of-freedom). In 
this article, we also go beyond the Born-Oppcnheimer ap- 
proximation to find the non-adiabatic corrections to the 
Berry phase. However we are interested in a two-level 
system (qubit) which is very weakly coupled to each of 
an enormous number of environment modes (or classical 
noise modes), such that the combined effect of all these 
modes on the system dynamics is very strong. This is 
the opposite limit from that considered in the works on 
molecular dynamics. 



II. BERRY PHASE DUE TO A 
NON-MARKOVIAN ENVIRONMENT. 

Here we avoid usin g t he adiabatic renormalization 
method of Leggett et al [13| ■ Instead we perform an exact 
"polaron" transformation on the Hamiltonian in Eq. l[8|). 
which fits with the spirit of adiabatic renormalization, 
while being a much more controlled approximation 14 1 . 
This elegant approach is standard for polarons [34| . and 
was first applied to the spin-boson model in a num- 
ber of papers [H, l3rl l37l . l38l l39ll . some of which arc 
much neglected. Refs. |37l l38t I391 showed that the non- 
interacting blip approximation [1 31 ] is given by a simple 
weak-coupling analysis of the transformed Hamiltonian. 
Elsewhere, we will present a detailed review of this ap- 
proach and discuss its regime of validity (greatly over- 
estimated in Ref. [39|); here we simply note the remark- 
able conclusion that the polaron transformation can map 
a spin coupled to a highly non-Markovian environment 
onto a spin coupled to a almost Markovian environment. 
The latter can then be treated with a Bloch-Redfield mas- 
ter equation. 



A. Transforming to the rotating basis. 

To deal with a problem in which the axis the envi- 
ronment couples to is rotating with time (as sketched in 
Fig- H]), we go to a rotating basis Q whose 2-axis re- 
mains parallel to e(t). We transform to such a basis us- 
ing U = ex-p(—i^ipa z ) exp(i^9a y ) exp(i^(pa z ), in terms 
of polar coordinates e(t) = (6, ip). This choice of U gives 
U(e(t) ■ a)W = <j z while having no ambiguity for 6 = 0. 
In this basis, the Hamiltonian has an extra magnetic field 
equal to the basis' angular velocity [l^], u>, thus it is given 
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by 

Kytbenv = i(dW/d*)W t + Z^ys&cnv^ 
"2 



(8) 



~T~ + %) + Wenv({a], aj}) , 



spin 



where the angular velocity in this rotating basis is 
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-cpsmtf coscp — t 
-ip sin 6 sin 99 + 6 
— (p (1 — cos ( 



sin p 

cos y> 



(9) 



For convenience in what follows, we define u± as the 
magnitude of the component of lo that is perpendicular 
to the environment-coupling axis, then 



2 2 
IjJ 1 = ol 



<^ 2 sin 2 1 



(10) 



Thus to summarizing the situation, we have removed 
the time-dependence from the coupling to the environ- 
ment, by going to the rotating frame. The Hamiltonian 
in Eq. JS} is a biased spin-boson model (with the envi- 
ronment coupling to the spin's z-axis) and thus we can 
proceed to treat it via a polaron transformation in a sim- 
ilar manner to Refs. M, M, iH US [U . 



B. 



Physics of the polaron transformation to the 
basis of shifted environment modes 



The physics behind the polaron transformation is the 
idea that fast environment modes tend to adiabatically 
follow the slow system modes, so they can be thought 
of as a "cloud" around the system, renormalizing its dy- 
namics. So we go to the basis of this "cloud" , by writing 
^sys&E m terms of a new (orthonormal) set of basis-states 
in which environment modes are shifted due to the force 
induced by the spin being f or L A cartoon of this (for 
a single environment mode) is shown in Fig. [5] 

In this basis of shifted environment modes, any wave- 
function can be written as 

IV>) = [ U hrn 1 ,-,m N U}\m\,--- ,m ] N ) 

{™j) 

,mi ,rriN 

|i)|mjy • • ,m\f) , (11) 



where indicates that the jth environment mode is in 
the mth eigenstate of the shifted Hamiltonian Tl^ = 
Kj(a]: + a,j)/2 + ilja-cij (we drop an irrelevant constant 



term). Similarly |mj 

>shift;| 



is the mth eigenstate of fi 



shift a 



given by W • ' with Kj 



—Kn. We assume that each 



environment mode is only weakly affected by the system, 
Kj -C Oj, then the difference between |mj) and | 
small. Dehning ctj = (Kj/flj), one finds, 



is 



(mj|mt) ~ l-(m 3 + l/2)c 



(12) 




Adiabatic 
spin-flip 

c=0 




env. mode j 



env. mode j 




Non-adiabatic 
spin-flip 



env. mode j 



FIG. 2: A cartoon of adiabatic and non-adiabatic spin-flip 
transitions with shifted environment basis-states. Initially 
the spin is aligned along the direction of the coupling axis 
(indicated here by z), the environment modes are all shifted 
slightly to the right by the interaction with the spin; here we 
only show the jth mode which is in the shifted ground-state la- 
belled \A) (the dashed-line indicates the equivalent unshifted 
state). We then show two possible spin-flip processes. The 
adiabatic spin-flip, involves the jth mode going to state \B), 
which is the ground-state of the environment mode shifted 
slightly to the left by the interaction with the spin. The non- 
adiabatic spin-flip involves the jth mode going to state \B*), 
which is the excited-state of the environment mode shifted 
slightly to the left by the interaction with the spin. 



where aj = Kj/flj. The higher the level, 



Tfij , of the 



initial state of the jth environment mode, the smaller 
this overlap is. While this overlap is very close to one, 
the fact that it is slightly less than one can have a huge 
effect on the system, because there are many such over- 
laps. Each time the spin- flips it must carry the "cloud" 
of environment modes with it, as a result the spin-flip 
rate is multiplied by the product over all j of the above 
overlaps; this product of N overlaps — each of which is 
slightly less than one — will decay exponentially with TV. 
Thus for large N, the A^-particle environment state after 
the spin-flip will be almost orthogonal to the TV-particle 
environment state before the spin-flip. This means that 
the matrix element for spin-flips is suppressed (becoming 
exponentially small in TV). Such a suppression is often 
called the Anderson orthogonality catastrophe [13] , since 
he pointed out that it can strongly suppress the tun- 
nelling of an electron between two metals electrode [4(| ■ 

Intriguingly the higher the level, rrij, of the initial state 
of the jth environment mode, the further the overlap in 
Eq. (JTSJ) is from one. This would imply that the orthog- 
onality catastrophe is stronger at higher temperatures. 
However this view is too simplistic, because it neglects 
possible excitations of the "cloud" of environment modes 
that follow the spin. The overlaps which correspond to 
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such excitations are 

((m-l)j|mj) ~ m) ,2 aj, 

((m + l)j|mj) ~ -(rrij + lY^aj. (13) 

Note that if we flip all the spins (T^l) in Eq. (fll5j) . then 
the right-hand-side changes sign. We will see that these 
effects tend to counteract the orthogonality catatrophy 
which occurs in the adiabatic evolution. Thus we must 
take seriously their contributions to the dynamics. 

C. Transforming to the basis of shifted 
environment-modes 

We transform Eq. §E§ to the basis of shifted 
environment-modes, using a polaron transformation [34l . 
l35l l3l| ; for completeness we explain the transformation 
in Appendix [Bl The transformed Hamiltonian is 

/T/shift nil I nil I -\ 1 

sys&cnv ^sys ' ^"env ' ' 

7i' sys = -\{u) z a' z +Lu x c~ F <r' x ), 

3 

V = -±u x (a' x Q even + fr' y Q odd ), (14) 

The non-interacting Hamiltonian, (H' sys +H' cnv ), involves 
no transitions of the environment modes; thus it gives 
adiabatic evolution. The interaction term, V, contains 
all non-adiabatic effects (i.e. transitions between environ- 
ment states). The first term in V involves transitions of 
an even (non-zero) number of environment modes, while 
the second involves transitions of an odd number of envi- 
ronment modes. The Hermitian operators Qcvcn.odd are 

Qcvcn = \e- F {U- + 11+), 
Qodd = ±c- F (K--K+), (15) 

where 7Z± = FJ ■ [l ± ajd-' z (a'j — a^)]. The exponential 
suppression of perpendicular fields in Ti' 8ys is given by 
the Franck-Condon factor, 

F=& = k^«i + ( 16 ) 

3 

If the environment consists of harmonic oscillators at 
temperature T, then F = ^ . a| coth (/?fi,-/2), where 
(3 = (hsT)- 1 (see Appendix EJ) . 

D. Berry phase from a Born-Oppenheimer 
approximation 

We can make a Born-Oppenheimer (BO) approxima- 
tion of Eq. (TT4")) . by simply neglecting V (since it ex- 
cites environment modes). Such a BO decoupling of 



fast modes (environment) from slow modes (system) is 
known to create Berry phases @. A priori, one might 
guess that this BO approximation is justified when the 
rotation-rate, u>, is small enough that lo <C Below 
we find the non-adiabatic corrections and show a poste- 
riori that in fact the BO approximation is valid when 
1/2 

lo <C il m Gdis ■ For the non-Markovian environments 
that we consider here (Gdis ^1), this condition is much 
less stringent that the above guess. 

In the BO approximation, the spin-dynamics are sim- 
ply given by Ti.' sys . Assuming the Franck-Condon factor is 
large, we can neglect all terms which go like e~ F , then the 
Hamiltonian is simply — jlo z &' z . The off-diagonal terms in 
this adiabatic Hamiltonian arc exponentially suppressed, 
due to the Anderson orthogonality catastrophe [l2| as 
outlined above. Thus the spin adiabatically follows e(t), 
and the total phase acquired by the spin is ±| f* p dtoo z 
where ± are for f, {. To observe this phase, we study the 
precession of a superposition of \ and J,, this precession 
is given by the phase difference between ] and |, which 
is ^totai = J P dtuj z This equals the Berry phase, 

$bp = id<p(l-cos6) = A, (17) 

where A is the solid-angle enclosed by the path of K(i). 
Following Berry, we can write this as ^adK, and use 
Stokes' theorem to give $bp = / b dS where b = Vk x a 
is Berry's monopole-field, (bK,O0,b v ) = {K~ 2 ,0,0). 

This BO analysis is sufficient to show that the spin ac- 
quires a Berry phase, and that there is no dynamic phase, 
<f>dyn- However to see the form of the non-adiabatic 
phases, we have to go beyond the BO approximation. 

III. MASTER EQUATION ANALYSIS OF 
TRANSFORMED HAMILTONIAN 

To go beyond the BO approximation, and include the 
effect of V, we use a master equation approach. The ex- 
act master equation for the evolution of the the system's 
density matrix in the shifted basis is 

d f* 

-p' sys = -i[H' sys ,p' sys (t)}- + I drEhp^-r)], 

(18) 

where S[r; p' sys (t — r)] is a super-operator acting on 
p' (t— t) and is given by the sum of all irreducible inter- 
actions between the system and environment, traced over 
all environment modes [HI, 112] (see also Ref. [HI). This 
master equation is clearly non-Markovian, and is equiva- 
lent to the Nakajima-Zwanzi g eq uation [43j. One gets the 
Bloch-Redfield equation fl9l |20| by treating all terms in 
the Eq. (|18p to second-order in the system-environment 
interaction, this involves replacing X by its Born approx- 
imation S Born , while treating p' sys (t — r) to zeroth order 
in the interaction. Since there are effectively two different 
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^PsysW = -i^sys.PsysC*)]- + / dr § cvcn [r; ^ ys (t)] 



drS odd [T;p' (t)], (19) 



where §even,oddb~; -X'sys] 



"iBorn 
J even, 



odd [K^[-r;^ y8 ]] with 
Xg ys being any system operator. The master equation 
now looks Markovian (it is local in time), however some 
(weak) memory effects are encoded 0, [4f| [46[ in S (see 
discussion in Ref. (HI). The super-operator Seven acting 
on an arbitrary system-operator X' is given by 

Seven [t] X S y S ] 

= 4B^([Ar n W + iA c Z"W] a> x (o)K(-T)x^ s 

+ [AZZ(r)~iA^(T)} X' sys a' x {-r)a' x {Q) 
-[A%™(t) -L4^(r)] a' x (0)KysK(-r) 
-[^W+i^r(r)] a' x {-r)X' sys a' x {Q) 



The equation for §odd[/?~; -^sys 
"even" — * "odd" and a' x 
gives the noise functions, ^even' odd 



(20) 



is given by Eq. (|2"0|) with 
a' y throughout. Appendix [Cl 
We assume the en- 
vironment is defined by a smooth Caldeira-Leggett J(f2)- 
function [l6l ]. with a typical characteristic frequency f2 m . 
This function is defined by 



dfij(^) 



(21) 



However it is convenient to actually work with a di- 
mcnsionlcss J-function; defined by j(x) = il^ 1 J(xQ m ). 
Then Appendix [C] shows that for a highly non-Markovian 
environment, upon neglecting terms that are of order 
exp[— 2F] smaller than the leading order, one has 



A%™(t) ~ A s Z(r) * iRe 



A%£(t) * A o as / d m (r) 
where we define 

Gdis 



, Im 



-G di3 [(n m r) 2 +i2xin m T] 



-G diH [(n m r) 2 +i2 X in m T] 



dx j(x) coth[/3f2 m a;/2], 



(22) 



(23) 



Xn 



environments (one coupling to a x and the other to a y ) in and 
V, we identify two such second-order terms; 5^°™ con- 
taining a pair of <2 C vcn-operators and S^™ containing a 
pair of Q dd-operators. Note that cross terms (with one 
Qeven- and one Q dd-operator) drop out when we trace 
over the environment because they contain an odd num- 
ber of environment raising/lowering operators. We write 
P ' sys (t-r) = c m *^ p' sys (t) e- in '^ = K^[-r;^ ys (i)], 
where this defines the super-operator K sys . Then the 
Bloch-Redfield master equation is 



Gdis 



dx j(x) x 71 



l + (-l) 



i 



-pn„ 



(24) 



The square bracket in x-n is simply one for all odd n, and 
is coth[/3f2 m x/2] for all even n. Thus by construction 
\2 = 1. We note that the Franck-Condon factor, F = 
XoGdis- 



IV. NON-ADIABATIC PHASES DUE TO A 
QUANTUM ENVIRONMENT 

The phase information is contained in the off-diagonal 
elements of the system's density matrix. We write this 
density matrix as p' sys = |(1 + s' z a' z ) + s' + a' + + s'_a'_ 
where <54 = (a' x =fia' y )/2 are the spin-raising and lowering 
operators. Physically s' z is the expectation value of the 
z'-axis spin-polarization, and s'± = (s' x ± is' y )/2 where 
s' x and s' y are the expectation values of the x'-axis and 
j/'-axis spin-polarization. 

We can rewrite Eq. (|19|) as a matrix equation for s'_^_, 
s'_ and s' z , by using the fact that s'^ = tr[<7^/9g ys ], for 



/' 



-,z. We multiply both sides of Eq. (19]) by a'J 



and take the trace. To do so we evaluate the following 
traces 



tr sys [ a'Ja x o' v o' x (-t)\ 
= -tx sys [a'^a' y a' v a' y {-T)] 



= trsy B [^^(-T)^ai] 
= -tr sys [a' f a' y {-T) a' v a' y 



(25) 



and 



tr sys [ a'l a' x a' x {-r) a' v 
- tr S ys[<3-^ a' y o-y(-r) a' v 



= ^syslrfK o' x {-T)a' x ] 
= ^syslrfK a'y{-T)a' y ] 
SO 

OS* (26) 



2RcfSl 



where for compactness we define S = e~ luJzT . We do not 
need the matrix form of the traces in Eq. ([25]) . because 
the presence of the minus signs causes the terms to cancel 
amongst themselves, so they play no further role in our 
analysis. 

Now we note that the matrix equation for s' + , s'_ and 
s' z , is actually three uncoupled equations. Thus to see the 
phase information we need only analysis the equation for 
s',, which reads 



dt + 

where uj 

Eq. mm' 



dTAZZ(r)e- 



(27) 



is given in Eq. (fTfJ)) . and ^4^™(r) is given by 
The integral over r can be written as 



dr^ n (r)e- 



\ 2Gdis^r 



—0J Z 



2Gdis^n 



(28) 
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where we define 



i(b) 



dr exp 



5G dis ((fi m r) 2 +i2( X i+6)fi I1 



(29) 



The real part of 1(b) can found exactly by writing the 
integral as one from — oo to oo and completing the square, 
the result is exponentially suppressed for large GdisXii 



Rc[/(6)] 



r l/2 e -G dia ( X i+b) 2 /2 
fim(G dis /2) 1 /2 ■ 



(30) 



In contrast the imaginary part has no exponential sup- 
pression with GdisXi- For large GdisXi we can drop the 
quadratic term in the exponent, and find 



Im[7(6)] 



1 



^ m G d is(xi + b) 



(l + OKGcisXi)- 1 ]). (31) 



Now we use Eqs. (|30I31[) to expand Eq. (f2"8")) in powers 
of u> z . Then we cast the right hand side of Eq. (|27|) in 
terms of powers of w±_ z , so it reads 



dt + 



+0[7i 2 ] 



(32) 



where 



7i 



Q-^Tr/Gdis] 1 / 2 , 72 = (2fi m G di sXi)" 2 - (33) 



The w^-term is beyond Bloch-Redfield method applied 
in this work. We will discuss it in detail elsewhere, here 
we simply note that dimensional analysis shows that it 
goes like 7 2 (it is a third order contribution to Seven, odd 
and thus contains a double time-integral whose integrand 
is dominated by times of order 71 ) . It is trivial to solve 
the above equation for s' + , the solution is 

S ' + (t p ) = ex P [-i$total - Aotai] s+(0), (34) 

where the phase $totai and the dephasing factor Aotai 
are given by 



^total 
Aotai 



df 



u z - 72 u\u z + C[7i] u± 



Noting that w Xj2 go like t , we immediately get 
$dy„ = 0, 



$ B p 



At to. 



A, 



NA 

(2) 
NA 



0, 



(35) 
(36) 

(37a) 
(37b) 
(37c) 



dt ( 72 c4^+0[ 7l Vi), (37d) 



where A is the solid-angle enclosed by the coupling-axis 
(see Section lrTD)) . Higher-order non-adiabatic phases will 
appear, however just like the 0[cj^J term in the integrand 

(2) 

of ^j^^, most of them are beyond the Born approximation 
used to get the Bloch-Redfield equation, Eq. (fT§j) . 

We draw attention to the fact that $^ is zero. The 
term in (-i$ to tai - Aotai) which goes like t p is purely 
real; thus it contributes to the total dephasing factor, 
Aotai, but not the total phase, <& to tai- 

Now we turn to the total dephasing factor Aotai, we 
split it in to dynamic, geometric, and non-adiabatic con- 
tributions (just as we did for ^totai)- We define dynamic 
dephasing as terms in Aotai that go linearly with time 
(thus the usual exp[— tjTi\ dephasing of two-level systems 
by weak-Markovian noise is dynamic dephasing) . We de- 
fine geometric dephasing (a term coined in Ref. [26j ) as 
terms in -D total that are f p -independent, while we define 
non-adiabatic corrections to dephasing as terms that go 
like t p to some negative power. Then we immediately see 



that D 



dyn 



D 



BP 



D 



(2) 
NA 



0. while 



-^NA — 



dt 7ie 



(3) 

There will also be a finite -D^A term beyond the Bloch- 
Redfield analysis. 

We assumed from the beginning that the Franck- 
Condon factor, F = GdisXo is sufficiently large that we 
can neglect those exponentially small effects that go like 
c~ F . Thus for a finite temperature environment it is 
natural to say that GdisXi is also sufficiently large that 



3 -G di „x?/2 2 



(38) 



D 



(i) 

NA 



oc c~ GdisXl / 2 is very small. Although this is not 
the case at very high temperature, since then xi ~~ ► 0. 
However for any finite temperature, as one increases the 
coupling, Gdis, one will rapidly reach the situation where 



Ana is 80 small that dephasing is given by D 



(3) 

NA 



t~ 



(which we believe does not have a similarexponcntial sup- 
pression). Then we can effectively neglect all dephasing 
in the situation where t p is long enough for accurate ob- 
servation of the Berry phase. By this we mean that when 

(2) 

t p is long enough to ensure that <Pj^ 

(3) (2) 

will also have Aotai ~ Aa < $na 
be tiny. 

With the aid of the definition of 71 and Gdis, we find 



BP 



1. 



D£i < <f>£i < 1, SO Aotai will 



<-(^^ + °[<™ 2 »- 1 ])if' (39) 

where TV is the number of environment modes and K 
is the coupling to each one (so N K 2 gives the strength 
of the coupling to the environment). In the derivation 
we assumed strong dissipation, Gdis ~ NK 2 /Q 2 n ^> 1, 
thus the second term (which comes from the w^-term) 
dominates. Unfortunately this is the term whose exact 
form we do not have. None the less we see it is sufficient 
to make the NK 2 as large as possible to minimize $^1- 
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V. NON-ADIABATIC PHASES DUE TO A 
CLASSICAL NOISE-FIELD. 

Here we consider the situation sketched in Fig. [ft), 
and discussed in Section II CI As shown in Ref. jig , 
this problem is equivalent to the above quantum prob- 
lem for an environment of harmonic oscillators, with K 2 
replaced by l3Qj(A 2 ), in the infinite temperature limit 
(3 = (k^T) _1 = 0. If we assume there is a continuum of 
frequencies in the classical noise, it is natural to define 



dOJ cl (Q)(---)n 



(40) 



the factor of ir is to make the analogy with Eq. (f2~Tj) . 
Then J c \(£l) can be thought of as a measure of the noise 
power at frequency f2. Once again we define a dimen- 
sionlcss J-function; j c i(x) = J c i(xQ m ) . Then we can 
get results for classical noise from the results we have for 
a quantum environment; we simply need to replace j(x) 
in all the formulae in the preceding section of this work 
with /3f2 m x j c \(x) and then take the limit /3 — > 0. Then 
we find that 



and 
Even n : Xn 
Odd n : \n 



1 



da; jd(x), 



dx j c i(x) 



^dis JO x 

f3Q. m [°° Ax j c \(x) 



lim 

/3=0 



G 



dis Jo 



(41) 

(42) 
0. (43) 



In this situation, we see that Franck-Condon factor, F = 
GdisXOj remains finite as T — > oo. Thus if we have strong 
non-Mar kovian noise, i.e. large Gdis, we can have large 
F for which the BO analysis in section III Dl shows us 
that rotating classical noise will induce a Berry phase in 
the spin. However for such classical noise, GdisXi — * 0, 
so there will be no exponential suppression of Re[J(6)] 
given in Eq. (|30|) . In addition we cannot use the result 
for Im [/(&)] given in Eq. (f3"Tj) . because it is only for large 
GdisXi- Instead we note that for GdisXi ~~ * 0, one has 



lim Im[/(6)] 



26 



o[b 3 



(44) 



Hence 71 remains the same as in Eq. (|33[) . but now 

72 = pfi^Gds. 2 ]- 1 . (45) 

From this we find that the total phase, $totai, is that 
given by Eq. (f3"5|) . Thus all phases (the dynamic phase, 
Berry phase, and non-adiabatic phases) are basically the 
same for the classical noise as they were for the quantum 
environment. They arc given by Eq. (|37|) with the only 
difference being that now 72 is given by Eq. (1451) . 



Turning to the dephasing factor, -Dtotai> the fact that 
Re[I(b)] is given by Eq. ([50|) with GdisXi = means that 
in place of Eq. ([38]) we have 



D (1) - 

^NA — 



dt 71 ujI 



(46) 



(2) 

We still have Z?dyn = -Dbp = -Dna = Since now 
D^j^ has no exponential suppression, it will dominate 
the dephasing, which will be non-negligible. To observe 
the Berry phase in the presence of this dephasing, one 
requires that iJ^I < 1. Taking luj_ 



t p we have 



D 



(i) 

NA 



[/ °°dnj c ,(n; 



[VNA? t F 



(47) 



The right-hand side contains the integrated noise power, 

which is of order NA 2 for noise with TV modes each 

with an amplitude A. This means that the longer the 

evolution takes, the less dephasing there will be. Once 

£p ?i (NA 2 )^ 1 / 2 , we will be able to measure the Berry 

phase without worrying about dephasing. 

The leading non-adiabatic correction to the phase is 
(2) 

$na> it S ocs n k° 



(2) 

NA 



n 2 

NA 2 NA 2 



0[{NA 2 )-^^. 



(48) 



This is much the same as for the quantum environment 
(despite the different form for 72), once again it is domi- 
nated by the second term which comes from the o^-term. 

Thus to have small $^a one nee ds large NA 2 . 

These conditions are less strict than the equivalent one 
for a conventional Berry phase observation. Since the de- 
phasing is weak, the accuracy with which one can observe 

(2) 

the Berry phase, is determined by how small <i>^ is. The 
total phase will equal the Berry phase to within 0.1% , 
whenever $^a ~ 10~ 3< ^bp, this requires that the noise 
is strong enough that N 1 / 2 At v > 10 3/2 ~ 31. In this sit- 
uation Aotai ~ £>na ~ 10~ 3/2 , so exp[-A. tai] ~ 0.97, 
which means that dephasing is indeed extremely weak; it 
only suppresses the signal is only about 3%. 



VI. FINITE LEVEL-SPLITTING IN THE QUBIT 

Thus far, we have assumed that 7i S ys&env does not con- 
tains a 7i sys -term, where a 7i sys -term would be one that 
contains system-operators but not environment opera- 
tors. Physically this means that the qubit levels would 
be degenerate if we turned-off the environment coupling 
in Eq. ([3]) or turned-off the noise in Eq. © . However it is 
natural to ask what happens to the effects that we discuss 
in a system where this degeneracy is not perfect. Indeed, 
this may be important for observation in qubits, where 
one must tune a succession of external gate voltages and 
fields until the qubit is as close to degeneracy as possible. 
This means that we can assume that 7i sys is small, but 
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not exactly zero (one only has an exact degeneracy if it 
is due to some symmetry of the system) . 

Thus we now briefly revisit the calculations that we 
performed throughout this work, adding a small H sys - 
tcrm. This means the Hamiltonian is 

Hsys&env = T~tsys(t) +Hint(t) +H C nv({a],a>j}), 

H sys (t) = -i«7-B(f), (49) 



with Hint given by Eq. (|3b[) . If we now do exactly the 
same transformation to go the rotating frame as in Sec- 
tion III Al We see that the resulting Hamiltonian is the 
same as that given in Eq. ||5J) with the following substi- 
tutions 



lo x + B x 



co x 

Uz -> u z + B z , 



Uly + By , 



(50) 



where B x , B y , B z are the components of B(f) written 
in the rotating frame, so B z is the component of B(t) 
parallel to the environmental coupling axis (defined by 
e(f)). Thus following through the calculation, we find 
that Eqs. (|35l36p are replaced by 



4> 



total — 



(If 



uj z + B z 



-72 [{u x + B x f + {toy + By) 2 } [u„ + B z ] 



+0[ 7l 2 ] [(«* + B x f + (uj y + By) 2 ] 3/2 ' , (51) 



D 



total 



7iC -G dlaX ?/2 ^ + Bx) 2 + {uJy + 



-<D[(u + B) 4 



(52) 



Now noting that u> x ,y,z goes like fp 1 , while B x ^ y z is in- 
dependent of f p , we can expand in powers of t~ . Then 
we find the phases 



dyn 



df B z , 



BP 



(1) 
NA 



(If 



u z + 0[~f 2 B 2 LU, 1 fB 2 Lj] 



0[ l2 B n lB] t~\ 



and dcphasing factors 



(53a) 

(53b) 

(53c) 
(53d) 



if one wishes to observe the Berry phase alone to high 
accuracy, one should ensure that 'I'dyn is sufficiently 
small; for example, for an accuracy of 1 in 10 3 we need 
$dyn < 10~ 3 . This provides a simple criteria for Bt p , 
which correspondingly say how close to degenerate the 
states should be in the absence of the environmental cou- 
pling or classical noise The requirement that Bt p <C 1 
corresponds to B <C w, thus the higher the power of B 
in Eqs. (|53l [54]) . the smaller the contribution. 



Turning to dcphasing, we see that there are now small 
dynamic dephasing and geometric dephasing terms. For 
finite-temperature quantum environments, these are all 
exponentially suppressed for strong environment cou- 
pling, but for classical noise (xi = 0) they arc only 
powerlaw suppressed. However one should not forget 
that we are interested in the regime where B <^ tv and 
hence £>dyn *C -Dbp <C £>na- Thus the dephasing will be 
dominated by m situations where .D^A ^ s n0 ^ ex P°" 

(3) 

nentially small, and be dominated by in situations 

where -D^A ^ s exponentially small. We note that there is 
finite geometric dcphasing here, therefore we should add 
this system (one coupled to strong non-Markovian noise) 
to the list of situations in which geometric dephasing oc- 
curs (non-Hermitian Hamiltonians [47| , weak Markovian 
noise [2(|, or ultra-slow noise [HI). However in each of 
these situations, the form of the geometric dephasing is 
rather different. 



The presence of finite B gives a finite $dyn, thus one 
might ask if it is still easier to observe this Berry phase 
than a conventional Berry phase. The answer is yes. 
In the conventional Berry phase the same energy pa- 
rameter, E, that controls the smallncss of non-adiabatic 
phase oc (-Ef p ) -1 , enters the dynamic phase cx Et p . Thus 
choosing E to minimize the former, maximizes the lat- 
ter. Here non-adiabatic corrections are controlled by one 
parameter (the strength of coupling to the environment, 
Gdis), while the dynamic phase is controlled by another 
(the magnitude of the field, B). Thus both these un- 
wanted phases can be independently minimized (by tak- 
ing large Gdis and small B). Indeed, we can relax the 
condition on the smallness of B, by reducing f p . This 
can be done without enhancing the non-adiabatic cor- 
rections by increasing the environment coupling (or the 
classical noise power). 



D 



dyn 



D 



BP 



df 7l C- G ^/ 2 [B 2 x +B 2 y], 



(54a) 



df 2 11 e- G *»xV 2 {B x u x + Bytuy), (54b) 



dt 7ie -^x?/2 (w 2 + 2) _ 



Our first observation is that the system develops a dy- 
namic phase, 'I'dyn which is proportional to Bt p . Thus 



Finally, we note that Eq. (|53cj) indicate that there is 
now a small leading-order non-adiabatic correction to the 
phase, $na- However this can be neglected under the 
above condition that B <g; uj, because it is much smaller 

(2) 

than fjj|, that we treated earlier in this article. There 
(54c) is also a even smaller modification of the Berry phase, 
however since it is so much smaller than $na ^ can a ^ so 
be neglected (indeed it would be almost impossible to 
observe). 
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VII. TOWARDS EXPERIMENTAL 
OBSERVATION 

The effects that we discuss here should be observable 
in experimental set-ups similar to those used to observe 
conventional Berry phases, such as the superconducting 
qubits experiments of Ref. [j| or the cold-neutron exper- 
iments of Ref. @ . 

The effects induced by strong non-Markovian classical 
noise should be relatively straightforward to observe. If 
one is to able to drive the three-components of the spin 
with arbitrary signals, then one can drive them with a 
man-made noise signal. One could follow Ref. @, and 
construct this noise-signal on a computer, thereby con- 
trolling all aspects of the noise. Importantly one can then 
choose the components of the noise to be correlated in the 
manner that indicates that the noise has a natural axis, 
which varies with time around a closed loop. This would 
be sufficient to show all the effects that we discuss asso- 
ciated with classical noise; noise-generated Berry phase, 
absence of the leading non-adiabatic correction, dephas- 
ing which goes to zero as we make the experiment time 
longer. 

To see the same effects for strong coupling to a non- 
Markovian quantum environment at low temperature will 
be much more challenging. At present, we cannot think 
of an example of an experimental system in which this 
could be performed. However a prc-requisite for a quan- 
tum computer is that one has a large number of qubits 
with controllable couplings. To observe the environment- 
induced Berry phase that we discuss, one probably only 
requires a few tens of qubits (much less than required 
for quantum computing). One qubit would play the 
role of the system, and would be tuned to it degener- 
acy point. The other qubits would play the role of the 
environment, and would adjusted to have a variety of fi- 
nite excitation frequencies. Then the coupling between 
the system qubit and environment qubit would be varied 
in time as in Eq. With such a quantum environment, 
one can enter the regime where dephasing scales like t~ 3 
rather than t~ 1 (because the t~ ^term has an exponen- 
tially small prefactor). However apart from this addi- 
tional suppression of dephasing (which is already small 
for classical noise), this environment induces the same 
effects as classical noise. Thus, since classical noise is 
easier to implement experimentally, we see no significant 
reason to try to implement the quantum environment. 



VIII. CONCLUDING REMARK. 

This work — on the physics of the orthogonality catas- 
trophe and the associated strong renormalization of spin 
dynamics by a highly non-Markovian environment — ar- 
rives at two intriguing conclusions. The first is that it 
appears that adding coupling to an environment (of the 
right type) may make it easier to accurately measure 
Berry phases, when naively one would say an environ- 



ment always suppresses such phase information. The 
second is that even classical non-Markovian noise induces 
the strong renormalization of spin dynamics necessary to 
generate a Berry phase; even though such renormaliza- 
tion is based on the intrinsically quantum effect of van- 
ishing overlaps between environment states. This means 
that the classical limit of this quantum effect, must cor- 
respond to a classical effect. Unfortunately the nature of 
the calculation in this article (a formal mapping from the 
quantum problem to the classical one) does not help us 
understand this classical effect. In future work we plan 
to study the classical problem directly, and find a simple 
picture of the classical effect that plays a role analogous 
to the T — > oo limit of the orthogonality catastrophe. 

Finally, we recall that our results for a quantum envi- 
ronment are valid for arbitrary temperatures if the en- 
vironment is made up of harmonic oscillators. However, 
as discussed in Section [I B[ the low temperature limit of 
these results apply for an arbitrary environment (a bath 
of two-level systems, etc). 



IX. ACKNOWLEDGMENTS 

My thanks go to R. Leone and L. Levy for inspiring 
this work, to M. Berry and P. Bruno for insights into geo- 
metric phases, to S. Florens for insights into spin-bosons 
models and adiabatic renormalization, and to R. En- 
glman for insights that led me to the polaron literature. 
Finally, I thank T. Ziman for his general advice on RG 
and his ongoing support throughout this work. 



APPENDIX A: THE AHARONOV-ANANDAN 
PHASE 

There exists another geometric phase which has no 
non-adiabatic corrections; the Aharonov-Anandan (AA) 
phase 15j . Thus it is natural to ask why we are working 
so hard to create a Berry phase with suppressed non- 
adiabatic corrections, when the AA phase already exists. 
Our response to this question is that the nature of the 
Berry phase makes it much more easier to calculate in 
most experimental situations than the AA phase. Here 
we explain what we mean by this. 

The Berry phase is given by the solid-angle enclosed 
by the parameters of the Hamiltonian, which are directly 
controlled by the experimenter. Thus for any given time- 
dependence of those parameters the Berry phase can be 
calculated with simple trigonometry. 

In contrast, the AA phase is half the solid-angle en- 
closed by those quantum states that return to themselves 
[TEl ]. While the concept of the AA phase is certainly 
mathematically beautiful, let us consider for a moment 
the amount of work the experimenters have to do to cal- 
culate the AA phase for a given time-dependence of the 
parameters of Hamiltonian. They must first completely 
solve the time-dependent Schrodingcr equation for evolu- 
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tion from time zero to time t p to find those states which 
return to themselves. For an arbitrary time-dependence 
of the parameters of the Hamiltonian, this cannot be done 
analytically, so they must numerically simulate the quan- 
tum evolution to find these states. Having done this, they 
must then find the solid- angle enclosed by these states. 
This involves recording their behaviour for all time from 
zero to tp. For arbitrary time-dependence of the Hamil- 
tonian, this typically involves re-running the numerical 
simulation to get the time-dependence of those states. 
Only then can one use trigonometry to calculate the solid 
angle enclosed to extract the AA phase. 

Of course, there are certain simple cases in which the 
AA phase can be calculated analytically e.g., spin- 
half in a field which is static from time zero to time 
t' , then instantaneously changes to a new direction, and 
then remains static for times from t' to t p . However these 



cases are exceptional. 

Thus one could say that the noise-generated Berry 
phase, that we present here, fills a gap between the con- 
ventional Berry phase and the AA phase. It is as easy 
to calculate as a conventional Berry phase, while having 
much smaller non-adiabatic corrections. 



APPENDIX B: TRANFORMATION TO THE 
BASIS OF SHIFTED BATH-MODES 

We label the eigcnstatcs of a z with s = ±1 =t, I, then 
a x = Y^ s I — s) (s\. If there were only one bath mode (the 
jth mode), then we would insert a resolution of unity of 
the form 1 = 53 m |m|)(m|| to the right and another of 

the form 1 = J2 n \ n J a ){ n J S \ to ^ e gi ym g 



= El- S >< s l £[K S >K S K>KI + |(m-l)7'>((m-l)7'|m$>K| 



Mi) 



+ \{m + l)7')((m + l)7 s |mj)(mj| + • ■ 
1 - a|(a 3 -at + a)aj)/2 - aja' z(j) (a] - a 3 ) + 0[a 2 af, a? a?, a*] 



where we define aj = (m + l) 1 / 2 |(m + a'j = m 1 / 2 |(m — 1)|) and 



(Bl) 



(B2) 



Note that we assume that the environment is extremely 
large, and so the coupling to any given environment mode 
is tiny enough that we can keep only lowest (second) 



order terms in 



Since we trace over the environment at 



the end, terms of C[a 2 a] 2 ] or of 0[a 2 a 2 ] can be neglected, 
because they can only lead to contributions of C[a|]. 

There are N bath modes, so we insert N resolutions of 
unity on the left and another N on the right. We then 
define 



N 



K = El- S >< s l IlK s )K 



iY 



1 

1 



5>I-«>H nK*)KI = i 



i=i 



-i 



JV 



5>I*K«I IlK>< r 

s j=l 



i o 

-1 



> (B3) 



where the second equality (equating the operators to the 
usual Pauli matrices) is only valid if we are working in 
the basis of shifted bath-modes. Following the same logic 



as for one bath mode, one gets 

a x = a' x exp[-F] JJ [l aj a' z (a] - a 3 )] , (B4) 

3 

where we use the fact that ay is very small (while N 
is very large), to write (l — a 2 (a,ja^ + a^a,j)/2) ~ 
^exp [ - o^{flja\ + dtoj-)/2] = exp[-F}. Further, 

the assumption of very large N means that (F) re- 
mains unchanged during the evolution. Thus we re- 
place F by F = (F) given in Eq. (fH))) . For an environ- 
ment of harmonic oscillators at temperature T, (d]aj) = 
exp[- / 50 J -]/(l-exp[- ( 9%]) where (3 = (feBT)" 1 , then we 
get the result below Eq. (fTB]) . 

Eq. (|B4p is easier to handle if we extract the spin- 
operators from the product as follows. We note that 
all terms with even (odd) numbers of ays, have an even 
(odd) power of a' z . Thus as a' 2 = 1, terms with even 
numbers of a^s contain no spin-operator, while terms 
with odd numbers of a^s contain a single a' z operator. 
We separate terms with even and odd numbers of a^s, 
by noting that all even terms in Yiji^ ~ x j] are given by 



Yij [1 + x j\)- while all odd terms are given 
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by the difference of these two terms. Thus we get 

Ox = o' x eT F + a' x Q cvcn + a' y Q od d, (B5) 

fry = 6-' y eT F + (TyQeven- v'xQodd, (B6) 



o, = a. 



where Qcvcn.odd arc the Hcrmitian environment operators 
given in Eq. (|15p . Writing the "universe" Hamiltonian, 
^sys&env; m terms of these results, we arrive at Eq. (TH| . 



APPENDIX C: EVALUATING TRACES OVER 
THE ENVIRONMENT 

The Born approximation of £ in Eq. p8[) is most 
naturally written with the system and environment 



operators in the interaction picture, thus we define 



K,y(- T ) 



■ T and Qovon,odd(-T) 



-iH' 



T Qcvcn,oddC Qn " T ■ Then Qovon,odd(i) are given by 
Eq. (T5J with a J -> e in ^a'J and a'j -> e'^a'J . We 



(B7) define the symmetric and asymmetric noise functions 



^sym _ 



itr 

'even 2 onv 



/tasym _ _:1 +T 
-"-c ven 1 2 env 



pLnv[2eve„(0),Qeve.i(-T)] + 

Pe"w[Qcve„(0),Qeve„(-T)]_ , (CI) 



where [X,y]± = X y ± yx. We define Ag£ and A^/ d m 
in a similar manner with "even" — > "odd" . We note that 
^r„,od y d m are all defined to be real. To find A™ d m , 



we need 



tr cnv [^ nV Qeve„(0)Q OV e„(-T)] = e"^ [ftf* (0) + K- (0) - l) ( - T ) + (-T ) - l)] . 

where Ti±(— t) = | fl^Li (l T OijaJe~ injT ± aja^e 1 ^ 1 "). For diagonal p™ v (i.e. a thermal state), 



tr 



cnv;j 



env;i/i i -I ~'t\ It ~'t — IO.-t i -I iiliT\ 

p (1 + Ctjdj - CtjClj ) (1 - (t , C J + o ,(/ ,(' J J 

I 



2 iy (2) p -«V _ a 2 (i) e in^ 



1 ">/ ( ' 



(C2) 
(C3) 



where we define 



and 



(2) 



for the trace over all 2V environment modes we have 



a'jCLj). Thus 



tr cnv [pg nv Q eV e„(0)Qeven(-r)] = e" 2F ( cosh[2/(r)] - 1) 



(C4) 



'Ma 



where /(r) = | Ej=i 1 
and we recall that ctj = Kj/Clj. This definition of /(r) 
is such that /(0) equals the Franck-Condon factor, F, in 
Eq. (IT6]) . To get the last line of Eq. (|C4|) . we note that 



< I ctj I <C 1 and hence use the approxi- 



mation ]jf(l - /j) ^ II j ex Ph/i] = exp [ - £\ /,•] . 
evaluate other traces in the same manner, we need 



We 



tr en v[^ nV 2odd(0)Qodd(-r)] = e 

tre„v[po nV Qodd(0)Qeven(-T)] = 



-2F 



sinh[2/(r)], (C5) 
(C6) 



The second result comes from the fact that Q ven an d 
Qodd respectively cause the transition of an even and 
odd number of environment modes. Thus their product 
involves the transition of a odd number of modes, so the 
trace is zero. 

Using the definition of the dimensionless J-function 
below Eq. (flM)) , we write 



dx j(x) 



+ c 



iH m r x — {3 Q, m x 



(C7) 



and then approximate f(r) by its expansion in r, 



exp 



Gdis(xO - iXl^mT - 5(^mT") 2 

+0[x 3 (n m r) 3 })} 7 (C8) 



where Xn is given in Eq. (|24[) . with X2 = 1 by construc- 
tion. This approximation is natural for a non-Markovian 
environment, since it has Gdis 3> 1 (see Section fl A[) . since 
we can assume that when we integrate over r, we do so in 

1 /2 

a saddle-point approximation. So long as X3 "C Gdis , 
the quadratic term in the above exponents ensures that 
c-^ r ) decays to be much less than e^ - 1 long before the 
(f2 m r) 3 -term becomes relevant (which is why we were 
able to drop it from the exponent). For an environment 
without low-frequency modes X3 ~ 1; thus Gdis ^> 1, we 
find that A™ d 7 are given by Eq. 

For classical noise X3 = 0, so the conditions under 
which we can apply Eq. (|C8[) are a little different. Since 
condition that Gdis ^ xi is trivially fulfilled, saddle- 
point approximation is applicable under the condition 
that we can neglect the quartic term in the exponent, 



GdisX4(^ i] 



It can be neglected if Xi "C Gdis- For 



noise without low- frequency modes, X4 
dition is satisfied when Gdis ^ 1- 



1. so the con- 
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